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We present a model for describing nuclear matter at finite density based on quarks interacting 
with chiral fields, a and 7r and with vector mesons introduced as massive gauge fields. The chiral 
Lagrangian includes a logarithmic potential, associated with the breaking of scale invariance. We 
provide results for the soliton in vacuum and at finite density, using the Wigner-Seitz approximation. 
We show that the model can reach higher densities respect to the linear-er model and that the 
introduction of vector mesons allows to obtain saturation. This result was never obtained before in 
similar approaches. 
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I. INTRODUCTION 

The problem of studying nuclear matter with chiral La- 
grangians is not trivial; for instance models based on the 
linear cr-model fail to describe nuclear matter already at 
p ~ po because the normal solution in which chiral sym- 
metry is broken becomes unstable respect to the Lee- 
Wick phase. In Ref. Q Furnstahl and Serot conclude 
that the failure of the cr-model is due to the restrictions 
on the scalar field dynamics imposed by the Mexican hat 
potential. A possible solution to this problem is still to 
use a linear realization, but with a new potential, which 
includes terms not present in the Mexican hat potential. 
A guideline in building such a potential is scale invari- 
ance 

In QCD the invariance under dilatation is sponta- 
neously broken due to the presence of the parameter 
Aqcd coming from the renormalization process. For- 
mally, the non conservation of the dilatation current is 
strictly connected to a not vanishing gluon condensate: 
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(F a (x)F^{x)). 



(1) 



In the approach of Schechter fg}, and of Migdal and 
Shifman [7| a scalar field representing the gluon conden- 
sate is introduced and its dynamics is regulated by a po- 
tential chosen so that it reproduces (at Mean-Field level) 
the divergence of the scale current that in QCD is given 
by Eq. fl}. The potential of the dilaton field is therefore 
determined by the equation: 
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where the parameter e vac represents the vacuum energy. 
To take into account massless quarks a generalization was 
proposed in Ref. Q, so that also chiral fields contribute 
to the trace anomaly. In this way the single scalar field 
of Eq. (J2J is replaced by a set of scalar and pseudoscalar 
fields {a,ir,(p}. 

It has already been shown that an hadronic model 
based on this dynamics provides a good description of 
nuclear physics at densities about po and it describes the 



gradual restoration of chiral symmetry at higher densi- 
ties [|| . In the same work the authors have shown a 
phase diagram, where the interplay between chiral and 
scale invariance restoration lead to a scenario similar to 
that proposed by McLerran and Pisarski in [Ioj |. This 
is not too surprising since the large N c limit explored in 
[Toj should be well represented in chiral models as the 
one discussed in Q. It is therefore tempting to explore 
the scenario presented in |9| at a more microscopic level. 

The new idea we develop in this work is to interpret 
the fermions as quarks, to build the hadrons as solitonic 
solutions of the fields equations as in JlA , \la] and, finally, 
to explore the properties of the soliton at finite density 
using the Wigner-Seitz approximation. 

Similar approaches to a finite density system have been 
investigated in the past [l3l - [l9| . A problem of those 
works is that the solitonic solutions are unstable and 
disappear already at moderate densities when e.g. the 
linear a- model is adopted 18j. We are therefore fac- 
ing an instability similar to the one discussed and solved 
when studying nuclear matter with hadronic chiral La- 
grangians. The first aim of our work is to check whether 
the inclusion of the logarithmic potential allows the soli- 
ton crystal to reach higher densities. The second and 
more important aim is to check whether the inclusion of 
vector mesons in the dynamics of the quarks can provide 
saturation for chiral matter. 

We should remark that the version of this model with- 
out vector mesons has already been studied at zero den- 
sity in [20|, but using a different technique to describe 
the single nucleon. Instead no calculation exists with the 
vector mesons. Moreover, both versions of the model are 
unexplored at finite density. 

The structure of the paper is as follows. In Sec. [TTI we 
describe the model we are using, in Sec. IIIII we review 
the Mean-Field approximation and the hedgehog ansatz 
while in Sec. IIVI we present the technique adopted for 
projection. Later in Sec. [V] we describe the Wigner-Seitz 
approximation used to mimic a system at finite density. 
Next in Sec. I VII and in Sec. IVIII we present our results, 
firstly for the soliton in vacuum both at Mean-Field level 
and by adopting a projection technique and then we show 
the results for the Wigner-Seitz lattice of solitons. Fi- 
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nally, in Sec. I Villi we present our conclusions and future 
outlooks. 



II. THE MODEL 

In a chiral quark-soliton model quarks are coupled to 
mesons in a chirally invariant way. Following Refs. 0- 
5, 9] we consider the Lagrangian: 

Co =^{i^d p - g„(a + in ■ t^ 5 ))iJj 

-V(<j>,a,w). (3) 



Here ip is the quark field, a and tt are the chiral helds and 
(j) is the dilaton field which, in the present calculation, is 
kept frozen at its vacuum value tfio. 

An extension of this model, already discussed in 
Refs. [2h5, ,9], is to add the dynamics of vector mesons 
and to incorporate the idea of universal coupling [2l| . 
This can be achieved by considering the vector mesons 
as massive gauge fields, following also the scheme of Ref. 
[22| . The new Lagrangian is given by: 



Cvm = 0(«7"9 M - gAo + i7V ■ T75) + dpl^^ ' ( p » + 75 " 4 '- 
- ^-fu^ + £{p»aD»(T + D p tt ■ D»n) 

+ \m 2 P (p» -P^ + A^- A") + -m* ci^w" 
-V{fo,w) (4) 



where lu^ is a vector-isoscalar coupled to baryon current, 
p p and A p are respectively a vector-isovector and an 
axial- vector- isovector fields coupled to isospin and axial- 
vector current. The covariant derivatives for the chiral 
fields and the field tensors for vector mesons read: 



The potential is given by: 
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where the logarithmic term generates from © . The first 
two terms of the potential are responsible for the breaking 
of scale invariance, while the second line is needed to 
ensure that in the vacuum (f> = <po, a = <7q and tt = 0, i.e 
it provides spontaneous symmetry breaking. The last line 
explicitly breaks the chiral invariance of the lagrangian, 
giving mass to the pion. 

In this work we assign the following values to the 
masses of bare fields: — 139 MeV, m p = tua = 776 
MeV and m u = 782 MeV. For the sigma field, since there 
are no experimental constraints, we use m a = 550 MeV 
and m a = 1200 MeV which are typical values used in 
nuclear physics if the sigma has to provide the interme- 
diate range attraction. We relate the quark-omega and 
the quark-rho couplings to the corresponding couplings 
with the nucleons. We keep fixed g p — 4 and we vary g u 
between 10 and 13. The pion-quark coupling constant 
g n will vary from 3.9 to 5. By varying the coupling con- 
stants in the ranges indicated above we will be able to 
tune the attractive and the repulsive interactions, in or- 
der to obtain saturation when studying the total energy 
of the soliton at finite density. The constants B and (f>o 
can be fixed by choosing a value for the mass of the glue- 
ball and for the vacuum energy e vac , while S — 4/33 is 
provided by the QCD beta function and it corresponds 
to the relative weight of the fermionic and of the gluonic 
degrees of freedom. Finally the constant Vq ensures that 
the potential energy is vanishing in the vacuum. 

As anticipated, we choose to keep the dilaton field 
frozen. This decision (which obviously simplifies the dy- 
namics of the system) can be justified by the results ob- 
tained in 0, Q , where it has been shown that at low tem- 
peratures the dilaton remains close to its vacuum value 
even at large densities. Therefore the degrees of freedom 
of our model are limited to quarks, chiral fields and vector 
mesons, similarly to the linear a-model whose solitonic 
solutions have been obtained in [ll|, [22| ■ The potential 
can be written in a simpler form and it reads: 



+ A p Ap u ). (5) V(a,ir) = \i(a 2 + tt 2 ) - \ 2 2 \n{a 2 + tt 2 ) - a ™> (8) 



The pion field mixes with the longitudinal component of 
the a\ [H, [24[ and so to ensure that the pion gets its 
physical mass once the two fields are decoupled, we need 
to introduce the constant: 
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Figure 1. Comparison between the logarithmic (solid line) 
and the Mexican hat potential (dashed line). 



Here the vacuum value <jq is fixed to be equal to f n = 93 
MeV. 

It is interesting to compare the logarithmic with the 
Mexican hat potential. In Fig. Q] it can be seen that 
in the case of the Mexican hat potential it is relatively 
easy to restore chiral symmetry by climbing the maxi- 
mum located at the center. This is not possible in the 
case of the logarithmic potential as long as the dilaton 
field remains frozen. Since only at large temperatures the 
dilaton field changes significantly @ we can expect that 
at large densities and moderate temperatures this model 
provides more stable solitonic solutions. This is a crucial 
question which will be investigated in our paper. 

An important point in our approach is that we aim at 
describing the dynamics of nuclear matter by incorporat- 
ing all the interactions already at a quark level. This is at 
variance with e.g. the approach of Ref . pol [26j where the 
vector field uj was introduced only at the nucleon level, 
but was not present in the dynamics of the quarks. 



III. MEAN-FIELD AND HEDGEHOG ANSATZ 

The Euler-Lagrangian equations that follow from the 
Lagrangian given in eq. @ read: 

[i$ - gn(cr + iir ■ T75) + g p - ■ (fi + ^4) ~ y</# = 0, 

dV 

(3d p D^a = -fapAp ■ tv - g^ifj - —, 

da 

dV 

pd^ir = /3g p (-p p x D"n + A^a) - ig^r^ - — , 

- d^P^v = g P v v + m 2 p p v , 

- d^A pv = g p a v + m 2 p A v , 
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(12) 



Here v v and a„ are the vector and the axial- vector cur- 
rents: 

t 

v v = p p x p^ + A p x A pv + /3tt x D v -n + ipj„-ip, 

(13) 



A p x p^ 



/3tt x D v a - f3aD u -K 



T 



(14) 



The fields equations for the model without vector 
mesons, provided by the Lagrangian in eq. ©, can be 
obtained by the previous ones by fixing g p = g u — 0. 

The previous equations are relations between quan- 
tum fields. The starting point of our calculation is the 
Mean-Field approximation, where mesons are described 
by time-independent, classical fields and where powers 
and products of these fields are replaced by powers and 
products of their expectation values. The quark spinor 
in the spin-isospin space is: 
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where the spinor Xh, defined as: 



Xh 



75 (K > 



satisfies the condition for the Grand Spin G = I 
G X h = 0. 



(15) 



(16) 

J: 

(17) 



The hedgehog baryon | B) — \N c q)\a)\ir)\w)\p)\A) is given 
by the product of quarks and coherent states of mesons 
fields and it corresponds to a linear combination of states 
with 1= 1/2 and J = 3/2: 



\B) 



E 

J MM 1 



(-Y +m CjSm,- Mi \J = I, M = -h, Mi). 



(18) 

An explicit expression for the coefficients Cj will be given 
in Sec. El 

It can be shown that the hedgehog state is one element 
of an infinite class of degenerate solutions of field equa- 
tions at Mean-Field level [27|. The profiles of the chiral 
and the vector mesons fields in the hedgehog state are 
given by the relations: 
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At Mean-Field level the mesons fields are classical and 
the differential equations governing their dynamics have 
to be supplemented by the appropriate boundary condi- 
tions. For the single nucleon case we impose the following 
boundary conditions to the fields: 



u'(0) = v(0) = 0, 
<(0) = h(0) = o, 
p(0)=w , (0) = ^(Q)=A T (0)=0, 

while at infinity the boundary conditions read: 



(19) 



a h {oo) = (t , h(oo) 
v(oo) 



0, 



-ga 



u(oo) V ~9 <7 o ~ e 
w'(oo) = p (oo) = A' s (oo) = A' T {oo) = 0, 



(20) 



where e is the quark eigenvalue. 

The total energy of the soliton at Mean-Field level is 
given by: 



Emfa = 47r 



Jr 2 dr(Ei, 



Ekin,Q + E a + E„ 



Euj + Ep + E A + E, 



pot) 



(21) 



where explicit expressions for each term in the energy 
density can be found in Appendix I. A test of the nu- 
merical accuracy of the solution originates from another 
way of expressing the energy, obtained by Rafelski [28[ 
by integrating out the fermionic fields: 
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(22) 



Our solutions satisfy this consistency test up to a preci- 
sion of the order of 10~ 3 . 



IV. PROJECTION 



il = (a, (3, 7) and Njmj is a normalization factor. 
Since the hedgehog states are eigenstates of the Grand 
Spin G, it is equivalent to rotate either in spin or in 
isospin space. Moreover, when studying diagonal matrix 
elements of nucleon states, it is customary to work with 
states where the third component of the angular momen- 
tum M is equal to —Mi since in this case the expression 
of the Wigner function is particularly simple. In this way 
the projection operator becomes: 



PjM = 



2J+ 1 



8tt 2 



d 3 nD J MM (n)*R(n). (24) 



The normalization factor has been determined by us- 
ing (E3J)-<I24J): 
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(25) 



Finally, the coefficients Cj in eq. (|T5)) are given by the 
expression: 



Cj = 



(B\P JM \B) 
2J+1 



8tt 2 



d 3 nD J M M (Q)(B\R(Q)\B). 



(26) 



Once we obtain the projected state, we proceed to eval- 
uate the corresponding energy. Basically we need to cal- 
culate the expectation value of the Hamiltonian on the 
projected state given by ([2"3"f. The projected energy can 
be written as: 

Ej ={JM - M\ : H : \ JM - M) 

=4tt J r 2 dr{E mt + E kintQ + Ej, a + Ej t „ 

+ E U + E.;. p + Ej. A + E J:POt ) (27) 

More details about the projection technique and the 
evaluation of the each term in the energy can be found 

in El HI. 



The hedgehog baryon defined in section IIIII is not an 
eigenstate of isospin and angular momentum. Moreover 
this semi-classical solution also breaks the translational 
symmetry of the Lagrangian, since the localized soliton 
is not an eigenstate of the linear momentum, either. In 
this work we restore the invariance under rotations by 
using the projection technique developed in [THHH. In- 
stead, we will not restore the translational invariance of 
our soliton. The spin-isospin eigenstates are defined as 
follows: 

\JMMi) = Nj Ml J d 3 nD^_ Ml {Q)*R(Q)\B) (23) 

where the weight functions D are the Wigner func- 
tions, i?(f2) is a spatial rotation through Euler angles 



A. Projected observables 

The formalism needed to compute most of the observ- 
ables can be found in [n], [29[. The only quantity for 
which we need to provide a new explicit expression is the 
potential energy, since in our potential there is a loga- 
rithmic term, not present in the a model. 
The matrix element for which we need to develop a full 
calculation is: 

E.j, po t = (JM - M\ : J d 3 rV(a h , h) : \JM - M) (28) 

More details about the calculation of this term are given 
in Appendix II. 
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To compute the static observables we have used the ex- 
plicit formulae given in the Appendix of Ref. |29] ; in ad- 
dition to radii and magnetic moments we also show the 
results for the average number of pions in the projected 
state, given by: 



(N^)j = N n C (J,N v ), 



(29) 



where Co (J, N„) [Tl| is a projection coefficient depending 
on the spin and on the average number of pions N n in the 
unprojected state. As already mentioned before, we do 
not perform a projection on the linear momentum, but 
we adopt an easier approach [3(| which provides a rough 
estimate of the center-of-mass corrections to the baryon 
total energy. The masses for the N (J — 1/2) and for 
the A ( J = 3/2) are then given by: 



Mj = (Ej - P 2 ) 1/2 • 



(30) 



V. 



WIGNER-SEITZ APPROXIMATION TO 
NUCLEAR MATTER 



In order to describe a soliton system at finite density 
we use the Wigner-Seitz approximation. This approach 
is very common and it has already been widely applied 
to both non- linear [TBI . [l9j and linear-a models . 
Specifically, the Wigner-Seitz approximation consists of 
replacing the cubic lattice by a spherical symmetric one 
where each soliton sits on a spherical cell of radius R 
with specific boundary conditions imposed on fields at 
the surface of the sphere. The configuration of the meson 
fields, centered at each lattice point, generates a periodic 
potential in which the quarks move. 

The spinor of quark fields must satisfy the Bloch the- 
orem: 



(31) 



where k is the crystal momentum (which for the ground 
state is equal to zero) and $fc(r) is a spinor that has the 
same periodicity of the lattice. 



A. Boundary conditions 

In the literature various sets of pos sible boundary con- 
ditions have been discussed [l8l. Qj|. In our work we 
adopt the choice of Ref. [l8[ which relates the bound- 
ary conditions at R to the parity operation, r — > — r. 
Respect to this symmetry the lower component v(r) of 
quark spinor, the pion h(r) and the rho p{r) are odd, and 
therefore they have to vanish at R: 



v(R) = h(R) = p(R) = 0. 



(32) 



Similarly, for the a field, the upper Dirac component, the 
lo and the A fields the argument based on parity provides 
the conditions: 

u'(R) = a' h (R) = uj'(R) = A' S (R) = A' T {R) = 0. (33) 
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Figure 2. Upper and lower components of the Dirac spinor 
((a), (b)), ct and pion fields ((c), (d)), in the model without 
vector mesons, as functions of the cell radius 7?. 
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- R = 2fm 

- R = 1 fm 

- R = 0.85fm 

- R = 0.74fm 




Figure 3. Vector mesons fields, as functions of the cell radius 
R. In (a) we plot the uj field, in (b) the p field, in (c) the As 
field and in (d) the At field. 



The boundary conditions at r = remain the ones 
given in eq. (|T9"]) . Basically the calculation consists 
in solving the set of coupled field equations in a self- 
consistent way for a given value R; practically we start 
from R = 4 fm, for which the periodic solutions are 
indistinguishable from the vacuum ones, and we slowly 
decrease the cell radius down to the smallest radius for 
which self-consistent solutions can be obtained. The pa- 
rameter set is the same used for vacuum calculations. 

In Fig. [3] we plot the Dirac and the chiral fields in the 
model without vector mesons for different values of R: 
down to R = 2 fm, the solutions do not change signifi- 
cantly, but as the cell radius shrinks to lower values, we 
see that all the fields are deeply modified by the finite 
density. 

For the model including vector mesons, we present in 
Fig. [3] the trend of the vector mesons fields. To better 
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— R=2 fin, Log.model with VM 

— R=1.2fm 

R= 2 fm, Log. model without VM 

-- R=1.2fm 




Figure 4. Color online. Baryon density as a function of the 
ratio r/R for the model without (dashed line) and with vector 
mesons (solid line). Two values of the cell radius 7? are shown, 
namely R = 2 fm and R = 1.2 fm. 



equations,which depend also on fc, as discussed in [18| . 
There it is shown that e top obtained by imposing the 
condition (|36|) is an upper limit to the top of the band 
and that the true top lies about half way between this 
upper limit and the bottom of the band. In conclusion, 
the band width turns out to be considerably smaller than 
the one estimated using the first method which leads to 
the upper limit etop in eq. (j34")) . 

Concerning the filling of the band, when working with 
chiral solitons at Mean-Field level the relevant quantum 
number is the grand-spin G and the lower band corre- 
sponds to G = 0. The only degeneracy remaining is color 
and therefore the three quarks per soliton are completely 
filling the band. The total energy of the cell is estimated 
by assuming a uniform filling and by averaging the energy 
within the band. 



VI. RESULTS: PROPERTIES OF THE 
NUCLEON IN VACUUM 



clarify the difference between the models without and 
with vector mesons, in Fig. 2] we show the baryon density 
profiles in the two cases. The relevant feature is that in 
the model without vector mesons the shape of the soliton 
becomes significantly more flat, at large densities, than 
in the case with vector mesons. This effect is due to 
the repulsion between the two solitons provided by the oj 
field, which prevents the baryon density to become large 
in the inter-nucleon region. This feature will have an 
important consequence on the dependence of the radii 
on the density, as discussed in Sec IVH Bl 

B. Band width 

We are solving a problem in which quarks are moving 
in a periodic potential and therefore Bloch theorem tells 
us that a band should form. How to define the width of 
the band is highly not trivial. The most sophisticated 
technique is the one provided in [l8j . 
Here we adopt two much simpler procedures. The first 
method is taken from [lj|, where the authors estimates 
the band width as: 

A = ^ +(^) 2 -bl, (34) 

etop — eo + A, (35) 

where eo is the eigenvalue of the ground state. 
An alternative approximation to the band width is ob- 
tained, following |lq |. by imposing that the lower Dirac 
component vanishes at the boundary: 

u(R) = 0. (36) 

A more accurate evaluation of the band effects can 
be obtained by solving in a self-consistent way the field 



We start by showing how the solitonic solutions for 
the single nucleon can be built and we check that it is 
possible to provide a reasonable description of the single 
nucleon properties with the chosen parameter set. In par- 
ticular we present here results for the lagrangian without 
vector mesons Cq for the set m a = 550 MeV and the typ- 
ical value g„ = 5 d, |Tl[ . Instead the parameters for the 
model with vector mesons are given by m a = 1200 MeV, 
g K = 3.9, g u = 12 and g p — 4. This set has been chosen 
in order to both get saturation at finite density and a rea- 
sonable description of the nucleon in vacuum. Anyway, 
since the calculation of nuclear matter properties done by 
using the Wigner-Seitz approximation is affected by large 
uncertainties we also present a second set of results with 
parameters better fitted to single nucleon properties. 

In Tables Q] and |TT] we present the static properties of 
the hedgehog baryon at Mean-Field level and we com par e 
them with results obtained in the linear a- model [22t l31j . 

In Tables ITTll and HVl we show the decomposition of the 
soliton total energy in its various contributions and again 
we compare with the linear u-model [22|, |3l| • We see that 
the contributions to the total energy coming from the 
chiral fields and from the vector mesons are comparable: 
as expected vector mesons play an important role in the 
dynamics of the soliton. It is also interesting to notice 
that the results obtained with the logarithmic model are 
very similar to the ones obtained with the Mexican hat 
potential, at zero density, specially when vector mesons 
are included. 

In Tables [V] and I VII we present the results after pro- 
jection in both models, with and without vector mesons. 
Moreover, in Table IVIII we present the results obtained 
using a second parameter set, better fitted to the single 
nucleon properties. It is important to remark that our 
results in the logarithmic model with only chiral fields 
are consistent with the ones obtained in 20]. There, a 
different approach based on the coherent pair approxi- 
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Table I. Various nucleon properties at Mean-Field level in the 
present work without vector mesons and in the cr-model [3ll | . 



Quantity 


LiOg. 1V10QG1 


(T-ivioaei 


r^xp. 


M(MeV) 


1176 


1136 


1085 


(r 2 )i=o 


(0.76 fmf 


(0.78 fmf 


(0.72 fm) 2 




3.83 


3.63 


4.70 


{rm)l=X 


(1.12 fmf 


(1.14 fmf 


(0.80 fm) 2 


9A 


1.27 


1.22 


1.26 




2.4 


1.9 


/ 



mation was used. Their results are similar to ours when 
the coherence length parameter x is taken to be of the 
order of one, as suggested in [32j . 

The results obtained both without and with the vector 
mesons in general overestimate the experimental values, 
particularly for the magnetic observables, once the pa- 
rameters are chosen so that the projected mass of the 
nucleon is close to its physical value. One has anyway 
to recall that for the mass an approximate correction for 
the spurious center of mass motion has been taken into 
account (see eg. ( [30]) ). while no center of mass correction 
has been done for the other observables. When this fur- 
ther corrections are taken into account the value of some 
observables typically reduces. 



Table III. Contributions to the soliton total energy at Mean- 
Field level in the Logarithmic model and in the Linear cr- 
model [31| . All quantities are in MeV. 



Quantity 


Log. Model 


Linear <r-Model 


Quark eigenvalue 


83.1 


107.4 


Quark kinetic energy 


1138.0 


1056.9 


E a (mass+kin.) 


OO4.0 




Ett (mass+kin.) 


486.0 


373.1 


Potential energy a — tt 


105.7 


120.7 


E qa 


-101.4 


-62.3 




-787.0 


-673.2 


Total energy 


1175.6 


1136.2 



A. Energy of the lattice 

In Figs. [5] and [6] we present the results of the total 
energy per unit cell in the present model and in the linear 
cr-model with and without vector mesons. For each given 
value of m a , the logarithmic model allows the system to 
reach higher densities. Notice that as m a increases, the 
system remains stable down to lower values of R because 
the chiral fields are more and more restricted to lay on the 
chiral circle. Moreover the introduction of vector mesons 



VII. RESULTS: SOLITONS AT FINITE 
DENSITY 

In this section we present the results obtained by 
studying a Wigner-Seitz lattice of solitons. We first dis- 
cuss the energy of the system at finite density and we 
then present the effect of the density on the single nu- 
cleon properties. 



Table II. Various nucleon properties at Mean-Field level in the 
present model and in the cr-model [22| with vector mesons. 



Table IV. Contributions to the soliton total energy at Mean- 
Field level in the Logarithmic model and in the Linear cr- 
model with vector mesons [22l ]. All quantities are in MeV. 



Quantity 


Log. Model 


cr-Model 


Exp. 


M(MeV) 


1329.5 


1331.7 


1085 


(r 2 e )i=o 


(0.78 fm) 2 


(0.76 fm) 2 


(0.72 fmf 


m=i (m^) 


4.49 


4.51 


4.70 


{rt)i=i 


(0.99 fm) 2 


(1.01 fmf 


(0.80 fmf 


9A 


1.34 


1.35 


1.26 


~Ntt 


2.62 


2.66 


/ 



Quantity 


Log. Model 


Linear cr-Model 


Quark eigenvalue 


114.5 


112.9 


Quark kinetic energy 


1075.8 


1080.6 


E a (mass+kin.) 


213.8 


212.2 


Ett (mass+kin.) 


393.2 


397.3 


Potential energy a — n 


81.2 


80.4 


E u (mass+kin.) 


-194.4 


-196.5 


E p (mass+kin.) 


162.6 


165.4 


Ea (mass+kin.) 


329.5 


334.1 


E qa 


6.54 


4.74 


Eqn 


-621.9 


-627.1 


Equj 


388.9 


393.0 


Eqp 


-163.8 


-165.9 


Eq A 


-341.8 


-346.4 


Total energy 


1329.5 


1331.7 
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Table V. Projected nucleoli properties in the present work 
without vector mesons and in the linear <r-model and com- 
parison with experimental values. 



Quantity 


Log. Model 


cr-Model 


Exp. 


E 1/a (MeV) 


1075 


1002 




M N (MeV) 


960 


894 


938 


E 3/2 (MeV) 


1140 


1075 




M A (MeV) 


1032 


975 


1232 


{r%) P (fm 2 ) 


0.55 


0.61 


0.74 


(r%)n (fm 2 ) 


-0.02 


-0.02 


-0.12 


(r 2 M )v (fm 2 ) 


0.7 


0.72 


0.74 


(rli)n (fm 2 ) 


0.72 


0.75 


0.77 


fip (jut) 


2.25 


2.27 


2.79 


/in (/Ltjv) 


-1.97 


-1.92 


-1.91 


9a 


1.52 


1.10 


1.26 




1.6 (J =1/2) 1.2 (J =1/2) 
2. (J = 3/2) 1.6 (J = 3/2) 


/ 



stabilises even more the solution, hence in comparison 
to the model with only a and n we can reach higher 
densities. 

In Fig. E] we plot the quark eigenvalue for the model 
without vector mesons as a function of the cell radius R. 
The line labeled e[^ p corresponds to the estimate of the 

top of the band given by eq. , while e[ b J p follows from 
eq. (|36l) . It is clear that in absence of vector mesons we 
never obtain saturation. Moreover to change the value 
of m, a does not modify this result. In Fig. [H] we show the 
analogous results in the case with vector mesons. The 
band structure is similar to the one in Fig. [SJ but here 
the main difference is given by a significant increase of 
the top of the band at high densities which allows us to 
obtain saturation. 
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Figure 5. Total energy of the soliton as a function of cell 
radius R for the linear cr-model plj and for the present model 
without vector mesons. Different values of m a are considered, 
= 5, g u = 12 and g p = 4. 
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Figure 6. Total energy of the soliton as a function of cell 
radius R for the linear cr-model [18| and for the present model 
with vector mesons. Different values of m CT are considered. 



To better understand our result which indicates the 
possibility of getting saturation, in Fig. [7] we compare 
the total energy of the soliton at Mean-Field without 
the contribution associated with the band. In order to 
emphasize the effect of the density on the energy we sub- 
tracted the mass of the nucleon in vacuum. It is clear 
that the exchange of vector mesons plays a crucial role, 
by contributing ^100 MeV at R — 1 fm, but it is not 
sufficient to get saturation. To determine which ingre- 
dients of the model are actually providing the repulsion 
at high densities, in Fig. \W\ we plot the interaction en- 
ergies for each term contributing to the total energy at 
finite density. The plotted quantities are defined as the 
value of the chosen energy contribution at R minus the 



Table VI. Projected nucleon properties in the present work 
and in the linear cr-model with vector mesons and comparison 
with experimental values. 



Quantity 


Log. Model 


cr-Model 


Exp. 


E 1/2 (MeV) 


892 


882 




M N (MeV) 


763 


750 


938 


E 3/2 (MeV) 


1030 


1029 




M A (MeV) 


918 


917 


1232 


(rl) v (fm 2 ) 


0.59 


0.58 


0.74 


(r|>„ (fm 2 ) 


-0.03 


-0.02 


-0.12 


{r 2 M ) P (fm 2 ) 


0.69 


0.69 


0.74 


(r 2 M )n (fm 2 ) 


0.70 


0.71 


0.77 


Hp (hn) 


2.72 


2.71 


2.79 


/in (/ijv) 


-2.49 


-2.5 


-1.91 




1.6 


1.48 


1.26 


(Nit) j 


1.1 (J = 1/2) 1.8 (J = 1/2) 
1.3 (J = 3/2) 2.2 (J = 3/2) 


/ 
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Table VII. Projected nucleon properties in the present work 
and in the linear cr-model with vector mesons and comparison 
with experimental values for the parameter set: g — 3.6, = 
13, g p = 4 and m CT = 1200 MeV. 



Quantity 


Log. Model 


cr-Model 


Exp. 


Ri in (MeV) 


1020 


1008 




M N (MeV) 


926 


912 


938 


£/ 3 / 2 [MeV ) 


1 1 A O 

114o 


11/17 
114/ 




]\/f * ( AJT 

Ma yivi ev ) 


lUOO 






(r%) P (/m 2 ) 


0.67 


0.66 


0.74 


<4)n(/m 2 ) 


-0.05 


-0.05 


-0.12 


(r 2 M ) P if™ 2 ) 


0.77 


0.76 


0.74 


{r 2 M )n (/m 2 ) 


0.78 


0.77 


0.77 


Hp (hn) 


2.63 


2.64 


2.79 


(UN) 


-2.37 


-2.38 


-1.91 


9a 


1.58 


1.46 


1.26 



corresponding vacuum value: 

Ei(R) = Ei(R) - Ei(oo). 



(37) 



From the figure it is clear that the band effect strongly 
influences the total energy of the soliton, by providing 
the largest contribution to repulsion at high densities. 
This is not surprising, because the band is associated 
with the sharing of quarks between nucleons. It is well 
known [33|, HH that in calculations of the N — N po- 
tential based on quark models the short-range repulsion 
is associated with the formation of a six-quark bag. In 



p -100 



-250 




Log. model with VM 

— ■ Log. model without VM 



R (fm) 



Figure 7. Total energy of the soliton at Mean-Field level in 
the logarithmic model without vector mesons (dashed line) 
and with vector mesons (solid line). The parameter values 
are: m a — 1200 MeV, g w = 5, g^ = 12 and g p = 4. 




Figure 8. Quark eigenvalue as a function of the cell radius 
R, in the model without vector mesons. The shaded areas 
represent the band as estimated in eq. (|34[) and in eq. (|36[1 . 
The first excited state 1 + is also shown. The quark mass in 
vacuum, here 465 MeV, is indicated by the dashed line. 



our calculation the exchange of vector mesons is the 
dominant effect at densities up to po, but at very high 
densities the band effect dominates. The total amount 
of these contributions leads to the repulsive mechanism 
responsible for the steep rising of the total energy at 
high densities, as it can be seen in the upper panel of 
Fig. m 

In Figs. [IT] and [T^] we show in details all the contri- 
butions to the total energy (without the band effect). It 
is clear that the uj meson provides short-range repulsion, 
partially compensated by the interaction of the quarks 
with the A meson. 

It has been discussed in the literature how to interpret 
the results obtained using the Wigner-Seitz lattice and 
in particular which should be the indications of quark 
deconfinement. In [l5| it has been suggested that decon- 
finement takes place when the upper band, which cor- 
responds to Grand Spin G = 1, merges with the lower 
band. In the case with vector mesons this occurs roughly 
at densities slightly larger than the saturation density. 
We should keep in mind that the estimate of the width 
of the band is affected by large uncertainties and it is 
well possible that in a more refined calculation satura- 
tion density and deconfinement density turn out to be 
well separated. 

Which conclusions can be drawn from our analysis? 
First and foremost we have a scenario in which it is pos- 
sible, for the first time, to obtain saturation by making 
use of the rather common idea of getting attraction from 
the chiral fields and repulsion from the vector fields. In 
this game the logarithmic potential plays a crucial role, 
by allowing the solitons to remain stable at densities large 
enough that the vector mesons can play a role. It is im- 
portant at this point to clarify whether we can get this 
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result only in a special and tiny parameters' range or 
if the mechanism leading to saturation is rather stable 
respect to the choice of the parameters' values. This 
analysis is shown in Fig. Q2] where we plot the value of 
the total energy, including the band effects, at the "sat- 
uration density" point minus the energy in vacuum for 
different values of the parameters, at fixed m a — 1200 
MeV. Here "saturation density" means the density at 
which a local minimum in the total energy appears, even 
if that minimum is not the global one. The minimum 
is global when the energy plotted in Fig. [T3] is negative 
and in that situation we are getting real saturation. In- 
stead, when the plotted energy is positive the local min- 
imum corresponds to a sort of metastable state. Since 
our calculation is affected by large uncertainties we think 
it is interesting to show also the parameters leading to 
this "false" minimum, because in a more sophisticated 
calculation (based for instance on a better estimate of 
the band) those energies could easily become negative. 
For parameters' values significantly outside the indicated 
range no local minimum exists. For instance for values 
of g u > 12 the local minimum of the energy disappears 
because the energy keeps raising as the density increases, 
while for small values of g u the repulsion cannot contrast 
the attraction and the energy gets smaller and smaller at 
large densities. It is important to notice that the range 
of parameters' values providing "saturation" is at least 
in part overlapping with the range of parameters' values 
for which a reasonable description of the single soliton in 
vacuum can be obtained. Finally, we have to recall that 
we are not really studying nuclear matter, but rather 
G = matter, which is composed of degenerate nucleons 
and deltas. Obviously, no experimental data exist for 



1400 




, 200 



m o =l200 MeV 
0=3.9, =l2,g =4 




Figure 9. Panel (a): total energy of the soliton with band 
effects (solid line) and without band effects (dashed line) as a 
function of the cell radius R for the model with vector mesons. 
Panel (b): the quark eigenvalue as a function of the cell ra- 
dius R for the model with vector mesons. The shaded areas 
represent the band as estimated in eq. (l34p and in eq. (|36|l . 
The quark mass, here 362.7 MeV, is indicated by the dashed 
line. 




Figure 10. Color online. Interaction energies Ei, as defined 
in text, as a function of the cell radius R in the model with 
vector mesons. Parameters as in Fig. [9] 



that type of theoretical matter, but we can expect it to 
be saturating, probably with a larger saturation energy 
than that of nuclear matter. The corresponding satura- 
tion density is also obviously unknown. 



B. Nucleon properties at finite density 

The question about the modification of the nucleon 
properties at finite density have been investigated in 
many anal ysis , both experimental [35] and theoreti- 
cal [26|, 130-39] . In Figs. Q3] and [15] we show how the 
values of a few observables evolve as a function of the 
spherical cell density p c = (4ILR 3 /3) -1 . The Mean- 
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Figure 11. Color online. Panel (a): contributions of the chiral 
fields to the total energy of the soliton as a function of the 
cell radius R in the model with vector mesons. Panel (b): 
contributions of the vector meson fields to the total energy of 
the soliton as a function of the cell radius R. Parameters as 
in Fig. I 
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Figure 12. Color online. Panel (a): interaction energies of 
the chiral fields with quarks as a function of the cell radius 
R in the model with vector mesons. Panel (b): interaction 
energies with the vector mesons as a function of the cell radius 
R. Parameters as in Fig. [9] 




Figure 13. Minimum of the total energy as a function of the 
density. Different sets of couplings are shown, for which the 
model admit saturation. See the text for more details. 



Field observables evaluated as a function of the density 
cell are the isoscalar electric and the isovector magnetic 
radius and the number of pions. The formulae for these 
quantities for the model without and with vector mesons 
read [11 El: 



r M/I=l 



r 4 (u 2 +v 2 )dr 



1 C 5 / s~~iQ /~~17T f~^P s~l A \ 

a*/=i j y 



(38) 



eld 



are given m 



and fii = i and the radial functions 
Appendix III. The problem with our evaluation is that in 
the Wigner-Seitz approach we are forcing a unit of baryon 
number in every cell of the lattice. Therefore, at densities 



large enough that the fields start occupying most of the 
cell and their value is no more strongly varying inside 
the cell, the various radii all simply scale with the size 
of the cell. The results we are obtaining are therefore 
indicative only at densities low enough that the fields 
are still relatively well contained inside each cell. As 
shown in Figs. [2] and [3] all the fields are well confined 
up to values of density close to po. This implies that the 
behaviour of the observables, evaluated from these fields, 
have no physical relevance at densities of the order or 
above nuclear matter saturation density pq. 

The problem we are facing is deeply associated with 
the Wigner-Seitz approximation in which the effect of 
the finite density to the various observables is due only 
to Hartree contributions. In the real case two neighbor- 
ing nucleons interact also via the Fock term. Notice that 
at densities large enough that the field fluctuations are 
suppressed, the Hartree contributions associated e.g with 
the pion field vanishes. On the contrary the Fock term 
becomes relevant at those densities and it provides contri- 
butions to the electromagnetic observables which in nu- 
clear physics are sometimes called pion-in- flight. These 
terms, of course, cannot be evaluated in the Wigner-Seitz 
approach. 

Another problem with the Wigner-Seitz approach is 
that it imposes spherical boundary conditions on the 
fields. This is particularly dangerous in the case of the 
chiral fields, since at Mean-Field level directions in ordi- 
nary space are connected with directions in isospin space, 
a situation which is certainly quite far from reality. Due 
to these problems a work is in progress 41| in which a real 
lattice will be studied, with boundary conditions which 
can change depending on the direction. We think it is 
nevertheless worthy to present our results in Figs. H~i]fl~5l 
so to compare them with future more precise estimates. 

For the model containing a and 7t only, our results 
can be meaningful at low densities where the dynamics 
is dominated by chiral fields. The introduction of vector 
mesons affects in quite a interesting way the isoscalar ra- 
dius: as shown in Fig. 2J the repulsion provided by the co 
field, for densities smaller than pr,, prevents the swelling 
of the nucleons. The qualitative effect of the inclusion 
of the vector mesons is to stabilize the shape of the soli- 
tons respect to compression. This can be seen also in 
the case of the magnetic radius where the reduction of 
this quantity as a function of the density is less marked 
than in the case without vector mesons. It is also pos- 
sible to evaluate the so-called "super-ratio" defined as 
(GE/G M ) p /(GE/G M ) vac , where G e ,m are the electric 
and magnetic form factors. In our calculation we obtain 
a reduction of this quantity as a function of the den- 
sity, similarly to what obtained in other works, although 
the effect here is much larger. Let us stress again that 
although the model with vector mesons allows to reach 
much larger densities, we do not attribute too much sig- 
nificance to the behaviour of the observables at p > p 
for the reasons explained above. In the lowest panel of 
Figs. [TJ] and Qj)] we also show the behaviour of the number 
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Figure 14. Isoscalar electric radius (panel (a)), isovector mag- 
netic radius (panel (b)) and average number of pions (panel 
(c)) as a function of cell density pc for the model without 
vector mesons. 
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Figure 15. Isoscalar electric radius (panel (a)), isovector mag- 
netic radius (panel (b)) and average number of pions (panel 
(c)) as a function of cell density pc for the model with vector 
mesons. 

of pions N-„ per unit cell at finite density. The possibility 
of an enhancement of the pion cloud, when the nucleon 
is not isolated, has been discussed in the literature [37j . 
In the present approach we obtain instead a decrease of 
Nv This result stems from the behaviour of the pionic 
field on the Wigner-Seitz lattice as shown in Fig. [5J The 
strong reduction of the pionic field is due to the bound- 
ary conditions requested by the spherical symmetry. It 
will be interesting to see if this result survives when a 
real lattice is used in the calculation. 



VIII. CONCLUSIONS 

We used a Lagrangian with quarks degrees of freedom 
based on chiral and scale invariance to study how the 
soliton behaves in vacuum and at finite density. We 
presented results for the simple model with just chiral 
fields and also for the model including vector mesons. 
To describe the single nucleon properties in vacuum we 
have used a projection technique. The values of the 
static observables in the logarithmic model at zero den- 



sity are comparable to the ones obtained using the linear 
er-model. 

For the description of the soliton at finite density we 
have employed the Wigner-Seitz approximation. We have 
shown that the new potential, which includes a logarith- 
mic term originating from the breaking of scale invari- 
ance, allows the system, for each given m a , to reach den- 
sities larger than the ones obtained with the er-model. 
Moreover, as expected, the addition of the vector mesons 
plays a double role in the study at finite density: on one 
hand it stabilizes the solution and allows to reach even 
higher densities, on the other hand it partially provides 
the repulsion necessary to obtain saturation. The re- 
maining repulsive contribution originates from the band 
effect. 

The possibility to obtain saturation seems to be a firm 
result of the model, at the level of the Wigner-Seitz ap- 
proximation. In fact, by exploring the space of param- 
eters, we have shown that the interplay between attrac- 
tion, provided by the chiral fields, and repulsion, given 
by the omega field and the band effect, allows to get satu- 
ration for a rather extended range of parameters' values. 

At sub-nuclear densities the dynamics should be dom- 
inated by the chiral fields and the modifications of the 
nucleon observables obtained in our work can therefore 
be physically relevant in the low-density range. In par- 
ticular the isoscalar radius presents a slight swelling, of 
the order of 5% in the model without vector mesons and 
at even smaller dependence on the density in the com- 
plete model. This trend is in agreement with previous 
calculations [M HI El ■ 

The present work will be extended in several directions. 
First a more precise and accurate calculation of the band 
in the soliton crystal will be done following Ref. [l8| . 
Work is in progress in order to study this same model 
by using the technique developed in Ref. [43[ , which pro- 
vides a more precise description of a multi-soliton system. 
Finally, the model can also be studied at finite temper- 
ature, including the dynamics of the dilaton field. We 
can expect that the effect of the finite temperature on 
the soliton lattice will be to reduce the stability, by low- 
ering the value of the dilaton field and therefore making 
it more easy for the chiral fields to fluctuate. It will be 
interesting to compare the obtained phase diagram with 
the one proposed by McLerran and Pisarski [10| . 
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APPENDIX I 

In this Appendix we provide the explicit expression for 
the energy density at Mean-Field given in eq. (j2"Tj) . 
The quark-mesons interaction and the quark kinetic en- 
ergies are: 

'* [g-nOhiu 2 - v 2 ) + 2g^huv + —uj(u 2 + v 2 ) 
^ 3 

2 



where: 



E-,, 



4tt 



-2g p puv + g p 



E. 



kin,Q 



3 a ( 2 1 2 

-As u z v l 

2 M 3 



dv du 



WA 1 



At: V dr dr 



(A.39) 



(A.40) 



and the energy density of the mesons fields and of the 
potential read: 



-^- 9p h(A s + 2 -A T 

dh (a 2 a 



h ( . 1 
h g p h + g p a h [ A s - - 



E u =--\ — 



E„ 



1 / duj \ 1 2 2 
, TO, ,CJ 

2 V dr J 2 w 



d P P 2 A 



A s - -A T 



(A.41) 



(A.42) 
(A.43) 

2 



1 

+ 2 



2 2 f„ 1 A 

-P- 9 P P ~9p \A S - -A T 



i 2 2 
+ TO„p 



pi 

(A.44) 



lcL4 7 



\ dr 3 dr 



— +g P P {A s + -A T 



-i 2 



\m 2 p (zAl + lA 1 



E 



pot 



--V((j) ,a hl h) 



(A.45) 
(A.46) 



APPENDIX II 



In this Appendix we provide a detailed calculation of 
the expectation value of the logarithmic potential given 
in eq. (|8]) between the projected states (see eq. ([27|)1. As 
already mentioned and shown in (Tl| , terms which do not 
involve the pion field (such as the quark-pion interaction 
energy and the a and quark kinetic energies) are not 
affected by projection. The main issue is the evaluation 
of the matrix elements of the chiral fields between rotated 
and unrotated hedgehog states. 
These matrix elements for the sigma field cr(r) are: 

(B\R{n)- x <r{r) n \B) = o{r) n {B^tiy^B) , (A.47) 



°(r) = \ (*(r) + £(Q)-V(r)) = a h {r) . (A.48) 

In an analogous way, for the pion field the matrix ele- 
ments become: 

(B| J R(^)- 1 7r(r)|B) = 7f(r)/i(r)(B| J R(fi)- 1 | J B) , (A.49) 



(A.50) 



Since the potential is a function of the pion only through 
quadratic terms, by using the previous relation we get: 

(Bii^nrV^rOiB) = g (n, e, MBiRiny^B) , 

g(n, 9, 0) = ±h(r) 2 (l + r- R^r) , (A.51) 

where the function g depends on Euler angles f2 and on 
the polar and azimuthal angles. 

For a generic function F of the quadratic pionic terms, 
the following relation holds: 

(BiRiny'Fi^irm) = Fig^e^miRiny'iB). 

(A.52) 

Therefore the projection of the potential term can be ob- 
tained by leaving the pure a terms unchanged and by 
replacing the quadratic terms of the pion with the func- 
tion g given in eq. (IA.51|) : 

A? K + h2 9(^ ^ 0) - A 2 ^ K + h2 9(M, 9, <j>)) 
- a a m\a h . (A. 53) 

The expectation value of the potential between the pro- 
jected states, eq. (|A.53|1 . becomes: 

Ej, po t = (JJ-J\ : J d 3 rV{a hl h,g{Q,e^)) : \J J — J) 



dr I suit 



1 

N J JO JO JO 

3 r>J in\Mt~ u „trt a / di S/oa — 1 1 



x J d&D^jityVia^h^i^e^MBlRiny^B) 

(A.54) 

where the Wigner function is equal to: 



2.7 



(A.55) 



Finally, the overlap between rotated and unrotated 
hedgehog states reads: 

(BiRinr^B) = 



p a + 7 

cos — cos 

2 2 



' -P^ + ^(4co S 2 fcos 2 ^-l N 



(A.56) 
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Here N a , are the average numbers of a and it mesons 
in the hedgehog state and Nj is a normalization inte- 
gral; explicit expressions of these quantities can be found 
in [ll|, [2l]. 



APPENDIX III 



The Mean-Field approximation for the isovector mag- 
netic moment reads: 



(A.57) 



where the radial functions G{ leld are given by: 
G Q M {r) J-uv , 



GJf(r) =4/3 



i/i 2 - g p ha h (a s - ^A T ) - g p h 2 p 



A s - -A T 



Gft(r) =(-1) 



4 I A' s - \A! T - ^f) ( As + \a t 



-Ag p p [A S + -A T 



(A.58) 



For the model without vector mesons, the coupling con- 
stant g p and the meson fields vanish. 
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